Electronic entanglement in late transition metal oxides 
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Here we present a study of the entanglement in the electronic structure of the late transition 
metal monoxides - MnO, FeO, CoO, and NiO - obtained by means of density-functional theory in 
the local density approximation combined with dynamical mean-field theory (LDA-I-DMFT). The 
impurity problem is solved through Exact Diagonalization (ED), which grants full access to the 
thermally mixed many-body ground state density operator. The quality of the electronic structure 
is affirmed through a direct comparison between the calculated electronic excitation spectrum and 
photoemission experiments. Our treatment allows for a quantitative investigation of the entangle- 
ment in the electronic structure. Two main sources of entanglement are explicitly resolved through 
the use of a fidelity based geometrical entanglement measure, and additional information is gained 
from a complementary entropic entanglement measure. We show that the interplay of crystal field 
effects and Coulomb interaction causes the entanglement in CoO to take a particularly intricate 
form. 

PACS numbers: 03.67.Mn,71.27.+a,71.15.-m,71.20.Be 
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Entanglement is a fundamental aspect of quantum me- 
chanics, responsible for a large range of complex phe- 
nomena not present in a classical setting. The entangle- 
ment of distinguishable particles was historically seen as 
something spooky, but is now considered a valuable re- 
source. It plays an essential role in quantum information 
theory, and has been studied in great detail both the- 
oretically and experimentally [H. Entanglement of indis- 
tinguishable particles has received less explicit attention, 
but it has been studied indirectly in both the quantum 
chemistry and condensed matter communities. Describ- 
ing the electronic structure of materials with a pure sep- 
arable state, represented as a single Slater determinant, 
is at the very heart of the modern computational ap- 
proaches, and the inability to do so is known under the 
term 'correlation'. The term encompasses both classical 
and quantum correlations (entanglement), where the for- 
mer results in mixed states and the latter in entangled 
states. The presence of entanglement is generally seen 
as a computational complication as it prevents the use 
of these standard approaches. The development of reli- 
able computational methods and entanglement measures 
are of key importance to turn also the entanglement of 
indistinguishable particles from a complication into a po- 
tential resource. 

Well known examples of strongly correlated materi- 
als arc the late transition metal monoxides (TMO) - 
MnO, FeO, CoO and NiO - which have been under in- 
tense experimental and theoretical attention for a long 
time[2(-y|. Here we report on a theoretical description of 
the late TMOs using the Local Density Approximation 
plus Dynamical Mean Field Theory (LDA-hDMFT) Q 
where the effective iinpurity model is solved by Exact 
Diagonalization (ED) [6| . The quality of the results is as- 
sessed through a direct comparison between the com- 
puted density of states (DOS) and photoemission exper- 
iments (XPS/BIS). We then take advantage of the direct 



access to the local many-body ground state of the impu- 
rity problem to analyse the entanglement in detail. 

The LDA-I-DMFT scheme is built around the mapping 
of the local lattice problem to an effective impurity prob- 
lem, which is then solved through various techniques. The 
ED solver is a natural extension of the Hubbard-I Ap- 
proximation (HIA)[8l-[l0|. In the latter the hybridization 
between the impurity and the self-consistent bath is ne- 
glected, leading to an atomic-like system described only 
in terms of the local projected LDA Hamiltonian H^^"^, 
a double counting term H^*^, and the on-site Coulomb 
interaction U. The Hamiltonian of the finite system is 
diagonalized numerically to produce an analytical self- 
energy. Such an approach gives excellent results for the 
spectral properties of actinidcs and lanthanides[^, [loj . 
However, in the TMOs the hybridization between the 
TM-3d and 0-2p orbitals is too strong to be neglected. 
While retaining the advantages of the HIA, the ED solver 
also takes into account a significant part of the hybridiza- 
tion. In ED additional auxiliary bath states are included 
in the local Hamiltonian, giving 
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where the indices i,j,k,l run over the local correlated 
orbitals and m runs over the auxiliary bath states. 
The energies and the hybridization strength Vim 
of the auxiliary bath states are determined from a fit- 
ting of the impurity hybridization function A(a;). The 
calculations were carried out in the paramagnetic (PM) 
phase, using a finite temperature (/3 = 0.00173 Ry) fully 
charge self-consistent LDA+DMFT implcmentation|10l- 
121 . Further technical details can be found in the Supple- 
mental MaterialllSl. 
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Figure 1. 
and NiO, 



(Colour on-line) PDOS of the TM 3d states (thick black lines) and O 2p states (dashed red lines) in MnO, FeO, CoO, 
and corresponding XPS/BIS data (black circles) 0]. The Fermi level is at zero energy. 



In Fig. [T] the calculated projected density of states 
(PDOS) of the TM-3d and 0-2p states are compared to 
experimental XPS (electron removal) and BIS (electron 
addition) dataQ showing mainly the contribution of the 
TM-3d states, due to the photon energies used. 

Both the occupied part of the MnO PDOS and the 
XPS data have a gradual onset developing into a broad 
structure between -1 and -8 eV, with the main peak at 
around -3.5 eV. The shoulder at -1 eV is well described 
by contributions from both TM 3d and O 2p states, and 
the sudden drop in intensity at about -8 eV is clearly 
reproduced. Notice that also minor experimental features 
like the satellite at -11 eV are found in our theory. The 
unoccupied part contains a single large peak at 5.5 eV, 
in excellent agreement with the BIS spectrum. 

The FeO 3d PDOS contains peaks at 1.5, -1.5, -4.0, 
-6.0, and -7.5 eV, all in acceptable agreement with the 
XPS spectrumQ. The peaks at 3.5 and -10 eV are not 
evident in the experimental data, which may be related 
to the fact that the experiment was performed with a 
non-stoichiometric sample (Feo.gsO), in contrast to the 
conventional FeO unit cell used in the calculation. 

The CoO XPS spectrum has a strong onset at low bind- 
ing energies, with a maximum at -1.5 eV, followed by a 



gradual decrease with a marked shoulder at -3.5 eV. The 
CoO 3d PDOS has two sharp main peaks at -1.5 and 
-3.5 cV, and an almost constant density from -4 to -10.5 
eV. The BIS spectrum shows a broad peak with a max- 
imum at 4 eV, which is located about 1 eV higher than 
the theoretical peak. The PDOS also contains a smaller 
peak at 5.5 eV, not visible in the BIS peak, likely due to 
the limited experimental resolution. 

The NiO 3d PDOS resembles the CoO 3d PDOS in 
that it has two main peaks at -2.0 and -3.5 eV, and a 
rather constant density from -5 to -10 eV, in excellent 
agreement with the XPS spectrum. It is even possible to 
assign the small features in the XPS spectrum at -5.5, 
-6.5, and -8.5 eV, to their respective peaks in the PDOS. 
The single peak at 4.0 eV in the PDOS is also in excellent 
agreement with the experimental data. 

Given the satisfactory comparison of the spectral prop- 
erties with experimental data, we now turn our attention 
to the entanglement in the thermal many-body ground 
state of the impurity problem, described by the density 
operator fF = e^^^ /Tr(e^'^^ ). There exist sev- 
eral ways to measure the entanglement in an N-electron 
many-body state[l[. In the case of a pure state \^), we 
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first look at the geometric entanglement measure [l^ 



i?G[|^)(*|] = l-max|(^''|vI/)|^ (2) 



where \^') is restricted to be pm-e and separable. From a 
computational point of view this measure is natural since 
l^"') corresponds to the best possible single Slater deter- 
minant description of the system. The search for the opti- 
mal 1 4'') is performed through a scan of the corresponding 
Hilbert space by recursively removing one electron at the 
time from the natural orbitalsfl^l- Although this proce- 
dure is in general not guaranteed to find the optimal I'f), 
it is exact for separable states and for pure 2-electron sys- 
tems, where it gives £'g[|*>(^'|] = 1 - Pmax with pmax 
the largest eigenvalue of the corresponding one-particle 
reduced density matrix pij = (\E'|c]ci|'I'). 

In order to analyse fP" it is necessary to generalize the 
entanglement measure Eq in Eq. ^ to mixed states. 
This can be achieved by replacing the overlap with the 
fidelity [iBl between and any separable state p', 



1 — max Tr 

p' 



(3) 



where p' = Y.^P^W^){'^'^l H^P^ = 1> and 1^-^) are sepa- 
rable states. Performing the restrained maximization in 
Eq. ^ is in general a formidable task, as one has to con- 
sider the effect of mixing several non-orthogonal separa- 
ble states. However, the problem can be simplified by not- 
ing that of a PM system with negligible spin-orbit 
interaction commutes with S^, and its ladder opera- 
tors S^. This implies that there is a common eigenbasis in 
which the many-body eigenstates l^*! „ ) can be indexed 
by s(s + 1) = (S^), and = (Sz), and that the eigen- 
values Ef are independent of m^. The thermal density 
matrix is therefore a mixture of several spin-degenerate 
eigenstates 



nis— — s 



where we have restricted ourselves to a single many- 
body spin quantum number s. We also observe that 
the spin-coherent operator Sa = exp(AS^)(l + |Ap)^^^ 
conserves the entanglement [l^ of the maximally spin- 
polarized state l^'f^s). Setting Am = exp(2i7rm/(2s + 1)) 
yields 



^ 2s + 1 



E 



2s + 1 



(5) 

where = 2Vy^(,.^| J (2s + 1). When applying the 

proposed Slater search algorithm to the pure state p^™ 
of the TMOs, the maximal overlap is always obtained for 



, where l^** /) is the closest separable state to 



l^fl^). Hence the separable state 
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P^ = 



2s + l 



(6) 



gives a local minimum of Eclpf]- Moreover each term 
p^™ contains the same amount of entanglement, origi- 
nating from the material specific l^'f^) and the action 
of the renormalization operator N^. Therefore we con- 
jecture that the local minimum given by p[ is in fact a 
global minimum. The geometric entanglement measure 
for the PM state p^ can then be written as [13 



I -Earn 
22''(2s + 1) 




(7) 



where p^ = X)iP»l*?,s)(^i,sl- Even when Eg[p^] is zero 
there is still a non-trivial term remaining in Eq. ([7]). This 
generic source of entanglement is very robust as it docs 
not depend on the details of the electronic structure of 
the system but only on (S^). 

As a complement to the geometric entanglement mea- 
sure we have used an entropic measure, based on the 
reduction of the many-body density matrix p to the one- 
particle reduced density matrix pij = Tr (pcjcj). This re- 
duction converts the entanglement in p to entropy, which 
can be measured directly [J e.g. in form of linear en- 
tropy Sl[p\ = Tr[p(l — p)]. However, quantifying entan- 
glement in form of entropy requires care, as any classical 
correlation in p is converted to entropy as well. We pro- 
pose the following entropic entanglement measure 

El [p] = Sl [p\ - min (Tr[p] , Tr[i - p])Sl [p] , (8) 

based on the fact that the entropy gained from the clas- 
sical correlations is less than the number of electrons or 
holes times the entropy in p[l3j. Removing this upper 
bound of the classical contribution simplifies the evalua- 
tion of El, but at the same time makes it less sensitive 
to detect entanglement in mixed states. 

Apart from the entanglement measures proposed 
above, other choices are possible. An intermediate scheme 
is to evaluate the geometric entanglement by using the 
relative von Neumann entropvfl^. However, the current 
lack of an efficient minimization procedure for the rela- 
tive entropy with respect to the separable reference state 
p' reduces its applicability to setting upper bounds on 



the entanglement [17 1 . 

The many-body eigenstate decomposition of p^ from 
the ED solver is shown in Table |I| The auxiliary bath 
orbitals are assigned a spin, but zero orbital angular mo- 
mentum. The ground state of NiO contains two sets of 
near degenerate states, with an energy difference of 80 
mcV. Although the energy difference is small, the low 
temperature strongly favours the eigenstates with the 
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Table I. Entanglement of the thermal ground state at T = 273K. The values of s is defined by s(s + 1) = (S^), and the 
degenerate eigenstate I'I'ms) corresponds to rUs = — s, .., s for each value of (Lz). E (eV) is the relative energy of I'i'ms), and 
Eg and El are defined in the text. 
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lowest energy. In all the studied materials, the ground 
state configurations maximize (S^) due to the strong 
Coulomb interaction. 

We start by looking at the entanglement in the (S^) 
and (L^) resolved components |'J'f„J('J'm, I of ■ As seen 
in TablelH both Eg and El increase in magnitude as |ms| 
becomes smaller. This trend is a direct consequence of the 
paramagnetism and is associated to the number of ways 
it is possible to distribute the electrons according to their 
spin. For such components Eg takes the form of 



EGimnJi^t 



1 



1 



EG[m){n])/{ 



2s 

s — nis 



(9) 



Let us now focus on the entanglement in the maxi- 
mally spin-polarized eigenstate l^'g). It depends mainly 
on the complicated interplay between the Coulomb in- 
teraction and the crystal field energies. The maximally 
spin-polarizcd ground state of NiO is dominated by a 
Ni high spin configuration (fiia where the subscript 
and the superscript stand for the number of and e^ 
electrons respectively). The on-site Coulomb interaction 
must preserve both (S^) and (L^) which implies that it 
can not couple this state to any other state. The Coulomb 
interaction is therefore effectively reduced it to a Hartree- 
Fock term, which makes the eigenstate IvPg) separable. 
The eigenstates at 80 meV have mainly tsil character. 
Here the Coulomb interaction is allowed to transform 
these states to fiii through pair-hopping. Nevertheless, 
when the maximally spin-polarized states are combined 
into the density matrix the entanglement in the states 
cancels out, giving Eg[p^] = 0.00. 

The ground state of CoO has mainly a Co tg^j con- 
figuration that couples to fiii through the pair-hopping 
induced by the Coulomb interaction. The pair-hopping 
gives rise to a strong entanglement in l^'^), and in con- 
trast to NiO, the entanglement is still present in , giv- 
ing Eg[p'^] = 0.02. 

The ground state of FeO has primarily a Fc con- 
figuration. The Coulomb interaction is again reduced to 
a Hartrcc-Fock term, except when the bath introduces 



an extra electron in the TM-3d orbitals. As a result 
the strength of the pair-hopping becomes very small, 
and it gives rise only to a very weak entanglement with 
i;G[|*'>(*'|] = 0.003 and = 0.012. 

The MnO ground state is dominated by the Mn fiio 
configuration. Even when one extra electron is introduced 
from the bath, the Coulomb interaction cannot induce 
any pair hopping, which makes p^ fully separable. 

The Coulomb interaction gives an additional contribu- 
tion to the entanglement when two electrons arc trans- 
ferred from the bath to the TM-3d orbitals due to 
the quadratic increase in repulsion energy. However, for 
TMOs this contribution is of the order of lO^"' for the 
geometric measure, i.e. too small to be seen in Table HI 

Let us now consider the entanglement in the thermal 
PM ground state p^. For NiO, FeO and MnO the Eg[p^] 
term in Eq. ([7]) is zero, and only the generic part con- 
tributes to the entanglement. A single Slater determi- 
nant method can therefore in principle obtain the max- 
imally spin-polarizcd states I'I'f j,), and then reconstruct 
the ground states through Eq. ([5]). However, such an ap- 
proach would not be adequate for CoO, as Eg[p^] is 
non-zero in this case. The generic contribution causes 
Eg[p^] to increase monotonically from NiO to MnO. A 
similar trend can be seen in the relative entropy between 
the locally projected [l^l thermally mixed density oper- 
ators in LDA and LDA±DMFt[i3]. However, due to a 
strong reduction of the large number of available many- 
body states in the LDA solution to only a few in the 
LDA+DMFT solution, the size of the relative entropy re- 
flects mainly the different degree of classical correlation 
in the two calculations. This trend in the relative entropy 
occurs in fact even when the LDA+DMFT ground state 
is replaced by the separable ground state of a correspond- 
ing Hartree-Fock-like LDA+U simulation. 

The proposed E^ measure was not able to resolve the 
entanglement in the PM phase, as shown by the nega- 
tive values in Table H] The condition used to subtract the 
classical contribution proved too strong for these systems, 
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and further refinement is needed to extend its applicabil- 
ity beyond pure states. 

Finally we note that the PM results can be extrap- 
olated to the type-II anti-ferromagnetic (AF) phase at 
zero Kelvin, by reducing the thermal ground state to the 
zero energy cigcnstatc with the largest magnetic moment. 
For NiO, FcO and MnO these eigenstates are separa- 
ble, which again makes them describable within theories 
working with a single Slater determinant, at least in prin- 
ciple. In CoO this state is rather entangled, with a fidelity 
entanglement of 0.12, which once more puts an upper 
bound on the accuracy of the single Slater determinant 
methods for this material. 

In conclusion, we have proposed a method able to sep- 
arate classical correlations from entanglement in a mate- 
rial specific theoretical framework. The ground states of 
the strongly correlated late TMOs have been studied in 
detail, and two main sources of entanglement have been 
identified. In the PM phase the entanglement is domi- 
nated by a generic contribution which is proportional to 
the atomic (S^). The second contribution, present also in 
the zero temperature AF phase, comes from the pair hop- 
ping induced by the on-site Coulomb interaction. This 
effect is supresscd by the presence of the crystal field 
splitting, but plays still an important role in CoO. We 
expect this behavior to be common to a wide range of 
materials with similar symmetries, e.g. Co doped ZnO. 
It would be of great interest to study the role of the en- 
tanglement in more exotic systems as complex actinides 
or Fe based superconductors, and see how it relates to 
their unconventional properties. 
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Supplemental material 

In the following formulas we have followed the conven- 
tion of using hats for many-body operators, e.g. U, and 
tildes for one-particle operators, e.g. A. A many-body 
operator B can be obtained from a one-particle operator 
B as 

B = ^B.,c^c,. (10) 

Exact Diagonalization 
Block diagonalization 

The Exact Diagonalization (ED) method solves the 
Anderson impurity problem by mapping it to a finite size 
problem, defined by the Hamiltonian 

H^^ = H" + U, (11) 

where H*' contains the one-electron terms, including the 
hybridization with a few auxiliary bath orbitals, and U 
describes the full Coulomb interaction between the elec- 
trons in the correlated orbitalsQ. 

An important technical issue in the implementation 
of ED is constructing the matrix representation of the 
Hamiltonian with respect to the many-body basis. This 
step is greatly simplified if the basis vectors are defined as 
single Slater determinants in some one-particle basis. The 
Slater determinants can be visualized in the occupation 
number formalism. For example 5 electrons in a 10 orbital 
manifold can form the following many-body states: 

\^\) = 11111100000), 
= 11111010000), 
1*^) = IllllOOlOOO), 

|^'5^_^) = 10000101111), 
= 10000011111). 

Here M is the number of possible many-body states, i.e. 
the binomial coefficient of 10 over 5. When the creation 
and annihilation operators are given in the same one- 
particle basis as ), then their action becomes a sim- 
ple remapping of the indices N — > N ±l,i — > j and a 
multiplication of a phase factor (±1). 

In a system with a large Hubbard U and bath ener- 
gies E far from the Fermi energy, it is usually sufficient 
to consider the ground-state configurations with N elec- 
trons, and the excited configurations with TV ± 1 elec- 
trons. The size of the many-body basis is then equal to 
all possible configurations of TV, and TV ± 1 electrons dis- 
tributed in K one-electron spin-orbitals. However, even 



with a moderate number of spin-orbitals {K = 30) and 
close to complete filling (TV = 0.8i^ = 24), the size of the 
Hilbert space becomes too large to handle in practice. 
However often the system under analysis possesses use- 
ful symmetries, which can be used to find some criteria 
to a priori identify a block structure in the Hamiltonian, 
and treat these blocks separately. 

A general way to obtain these criteria is to define the 
many-body basis vectors as the eigenvectors of a set of 
some commuting observables {A*^}, and then determine 
how the states mix under the action of the Hamiltonian. 
Since the basis vectors should be representable by sin- 
gle Slater determinants, only commuting one-electron ob- 
servables 

A'= = ^A^ctc„ (12) 

need to be considered. Furthermore the additional con- 
dition that the observables should commute with U al- 
lows to keep the calculation of the mixing rather simple. 
These two restrictions reduce the list of potential observ- 
ables {A*^} for the correlated orbitals to Sz and (or 
alternatively Sx and L^; and Sy and Ly). The auxiliary 
bath orbitals are not directly affected by U, so each bath 
spin-orbital m can be assigned an observable n"' = cJ^Cm 
that measures its occupation. 

Each spin-orbital j can be transformed into a common 
eigenstate of all these observables and assigned a vector 
of eigenvalues . A many-body basis vector de- 
fined as a single Slater determinant with respect to these 
orbitals, is trivially an eigenstate of the observables in 
{A''}, with an eigenvalue vector 

^^' = E(*f|6]£,l*f>a^'- (13) 

Obtaining these eigenvalue vectors is not enough to de- 
termine the block structure of , as the one-electron 
Hamiltonian H° does not in general commute with the 
observables in {A'^}. In particular, the hybridization 
with the bath orbitals rarely commute with . The off- 
diagonal elements of H° determine how the blocks will 
form. Each non-zero off-diagonal element H^„c}j^c„ al- 
lows a many-body basis vector to couple to a basis 
vector j) if the following condition is fulfilled 

P - A' ^ aJ'' - a'\ (14) 

The problem of generating a block structure can now be 
transformed into finding the connected components of 
an undirected graph, where the vertices are defined by 
{A-} and the edges by {a™ - a",0} through Eq. (fTi)) . 
This problem can be solved efficiently through the use of 
sparse logical square matrix multiplication, in the follow- 
ing steps: 
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1. Map each vertex to a matrix index, and each edge 
to a true element in the logical matrix T . 

2. Multiply T with itself, until it remains constant. 
This procedure makes the connected components 
complete. 

3. The true elements in a row or column in T gives the 
indices of all the vertices of a connected component. 

This algorithm can be improved by contracting dense 
regions of the graph before T is defined. Furthermore, 
the second step can be performed even more efficiently 
through the use of an update matrix V: 

1. V 

2. While V 

(a) T .or. V 

(b) V ^ (.not. T) .and. T.V 

Finally the Hamiltonian can is ready to be block diago- 
nalized by grouping all the many-body basis vectors with 
quantum number vectors {A^} matching a given con- 
nected component. 



Correlated orbitals and the hybridization function 

The ED method is built around fitting a few hybridiza- 
tion parameters V and E to the hybridization function 



wl-H 



LDA-DC 



(15) 



where G"(a;) is the bath green's function and '^^da-dc 
is the projected LDA Hamiltonian with the double count- 
ing removed. In our implementation the hybridization 
function is given on the Matsubara axis lo ~ iun = 
i-KT{2n -\- 1), and the fitting is performed by minimizing 
the cost function 



i^(V,E) 



vt[ 



a-E]-iV- A(c^) 



(16) 

through the conjugate gradient method. Here {W„} is a 
set of weights and ||A|j|, = Tr(A'''A) is the Frobenius 
norm. For the ED method to give physically relevant re- 
sults it is important that the original hybridization func- 
tion can indeed be approximated by a small number of 
peaks. 

The choice of correlated orbitals determines the map- 
ping from the lattice problem to the effective impurity 
model, and therefore the value of A(w). In order for the 
mapping to be realistic, the orbitals should be localized 
and nearly dispersionless. In the ED solver it is advanta- 
geous if the orbitals are eigenstates of the operator , as 
explained in previous section. Our current LDA+DMFT 



implementation [ll|, [l^ contains two choices of orbitals 
that meet these conditions to a large extent, the heads of 
the LMTO's (MT) and Lowdin's natural orbitals (ORT). 
The MT orbitals are given by 



(17) 



where is a step function, Yim is a spherical har- 
monic function, and 0i(r) is the solution to the radial 
Schr odinger equation within the mufhn-tin sphere of ra- 
dius S. The MT orbitals are both localized and eigen- 
states of \jz by construction. The ORT orbitals at a site 
R are defined as 



2. ^ 

keBZ j 



(18) 



where ^h^^^'^ir) is an LMTO orbital, O is the LMTO 

kj 

overlap matrix, the index j runs over all the LMTO or- 
bitals, and k runs over all the k-points in the Brillouin 
zone. The ORT orbitals are in general less localized and 
only approximately eigenstates to L^, due to the inclu- 
sion of the LMTO tail functions and the mixing through 
0~^l^ . However, when ipf^mji'^) correspond to localized 
states, both choices lead to similar orbitals. As an exam- 
ple, in the upper panel Fig.[3]the Ni-3d projected density 
of states of NiO in LDA is reported. The densities for the 
MT and ORT bases are practically identical. Neverthe- 
less, the small differences in the projections do have a 
large impact on A. In the lower panel of Fig. [2] we show 
the corresponding spectral density of the hybridization 
function 



(19) 



in ORT and MT. A striking difference between them is 
that the MT orbitals give rise to a huge high-energy peak 
structure starting at 10 eV, while this is nearly absent 
when the ORT orbitals are used. Such feature is common 
to all the transition metals oxides of the present study - 
NiO, CoO, FeO, MnO - and to many other systems as 
well. The strong high-energy hybridization given by the 
MT orbitals is very detrimental to the fitting process, 
since it attract the poles of the fitted function in Eq. (|16p . 
Therefore even though the ORT orbitals are less localized 
and require the use of a simplified basis [ll. 12 1, they are 
still to be preferred. The use of the ORT orbitals improves 
the situation, but it does not completely eliminate the 
high-energy hybridization. 

The weights Wn can be chosen to make the cost func- 
tion focus more on the physically relevant energy scale. 
In the calculations presented in the next section, Wn was 
set to sample a logarithmic mesh from 7 to 200 eV. The 
initial tight spacing of the mesh points gives more weight 
to the residues on the eV scale, while the sparsely spaced 
high-energy points make the cost function smoother and 
increase the stability of the procedure. The Matsubara 
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Energy (eV) 



Figure 2. (Colour online) Upper panel: The Ni-3d projected 
density of states of NiO in LDA. The black lines correspond 
to the MT orbitals and the red lines to the ORT orbitals de- 
fined in the text. Lower panel: Corresponding spectral density 
of the Ni-3d hybridization function in the same notation as 
above. 



frequencies smaller than 7 eV were excluded since they 
biased the fitting procedure in favour of a detailed de- 
scription close to the Fermi energy, at the expense of the 
overall description on the eV scale. 

The weights improve the fitting of the hybridization 
function on the Matsubara axis, but further refinement 
of the fitting procedure is necessary to obtain a good 
physical description for real energy values. To this end we 
used six bath states per spin-orbital in Eq. (fT6|). but only 
the physically relevant elements with the largest value of 
the utility function 



M(V,E) 



E 



-3|E|3 



-VV 



(20) 



were included in the final ED calculation. Here Af (V, E) 
is given in a basis where E is diagonal. The utility 
function M (V, E) excludes contributions from the broad 
high-energy structure described above, and improves the 
stability with respect to an initially metallic density. 

The post-selection procedure on the basis of the util- 
ity function worsens the fit on the Matsubara axis, but 



improves it closer to the real energy axis. It is important 
that a reasonable description of the original hybridization 
is obtained for both axes, and such situation is usually re- 
alized when the discarded states are well separated from 
the selected bath states. 



Computational details 

The TMOs were all set up in the NaCl crystal struc- 
ture. The lattice constants were taken from experimen- 
tal data, and the Brillouin zone was sampled through 
a conventional Monkhorst-Pack mesh of 12 x 12 x 12 k- 
points. The calculations were carried out in the paramag- 
netic phase, using a finite temperature fully charge self- 
consistent LDA-I-DMFT implementation[l3, [HI based 
on the full-potential linear muffin-tin orbitals (LMTO) 
method [l^. The double counting correction was defined 
as H-^*^ = fi^'^'N, where N gives the number of elec- 
trons in the correlated orbitals. /i^*^ was treated as a 
free parameter to ensure that the system has the cor- 
rect number of electrons in the ground-state and that the 
chemical potential is placed in the band gap according to 
the experimental photoemission spectra. 

The Coulomb interaction was parametrized by the 
Slater parameters F^, and F'^. The latter two were re- 
calculated at every iteration from radial integration of the 
bare Coulomb interaction and then multiplied with the 
screening factors 0.82 and 0.88, respectively. The screen- 
ing factors were set to reproduce the RPA screened Slater 
parameters for NiO in Ref. The parameter F^, i.e. the 
Hubbard U, could not be treated in the same way as it 
is too strongly affected by the screening, but was set to a 
fixed value from the start. The final self-consistent values 
of the Slater parameters are shown in Table iHl 

The ED calculations of NiO, CoO and FeO were per- 
formed with 10 TM 3d spin-orbitals and 20 auxiliary 
bath spin-orbitals. For MnO only 10 bath spin-orbitals 
were used, due to the greater computational effort. The 
final energy and hybridization strength parameters for 
the bath states are reported in Table [Hi Note that the 
values are given in the crystal field basis, e^ and t2g, in 
which the hybridization function is diagonal. The auxil- 
iary bath states between -17.9 and -19.5 eV correspond 
closely to the O 2s orbitals, while the rest have mainly O 
2p character. 



Entanglement 
Correlation and projection 

The term correlation has come to represent all elec- 
tronic or quasi-particle interactions that go beyond a 
simple single Slater determinant description. It can be 
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Table II. Self-consistent Slater parameters and auxiliary bath 
state fitting parameters. F" was held fixed at the tabulated 
value, while F^, F^, and the bath state parameters (Bath 
par.) were recalculated at each new iteration. The bath state 
parameters are given in the crystal field basis eg and t2g in 
which the hybridization function is diagonal. All values are 
given in eV. 



where 





Slater par. 


Bath par. (e^) 


Bath par. (t2g) 


TMO 


F« 






E V 


E V 


E V 


E V 


MnO 


6.0 


9.0 


6.1 


-5.1 2.1 




-6.2 1.4 




FeO 


6.5 


9.2 


6.2 


-4.4 1.9 


-17.9 2.1 


-8.9 0.7 


-4.6 1.2 


CoO 


7.0 


10.0 


6.7 


-4.6 1.8 


-17.9 2.0 


-7.5 0.9 


-4.0 1.0 


NiO 


7.5 


10.1 


6.7 


-5.7 1.9 


-19.5 2.0 


-8.8 0.9 


-5.1 0.8 



divided into two parts, classical and quantum correla- 
tion. The former gives rise to mixed states, and the lat- 
ter to entangled states. The four different combinations of 
classical and quantum correlation result in the following 
classifications of an N-electron state in terms of wave- 
functions and density operators 



N 



p' = 
p = 



pure separable > Uncorr. 



N 



Uij Y[ c]^ |0) pure entangled 
j p 

^p-|*'j)(\E'-| mixed separable 

i 

mixed entangled 



> Corr. 



where |0) is the vacuum state. Here pi and p[ are sets 
of probabilities normalized to one, and Uij are expansion 
coefficients. The vectorial indices i and j are composed of 
the ordered components ip and jp referring to a generic 
one-particle basis. For example the state cjcgCylO) cor- 
responds to a many body index i = (1, 3, 7). 

Classical correlations can derive from physical pro- 
cesses, e.g. thermal decoherence, but also be induced by 
a local projection upon some set of orbitals, e.g. the d- 
or f-orbitals of single atom in a solid. To see this, let us 
consider a set orbitals A labelled by the indices {jp}A, 
and define the operator subspaces 

^' = At = J2^^E^bjp e {jp}A,Y.\^A' - 

The full Hilbert space H can be partitioned as 



'H = A®B- 



A=[k^Q)-k^ e A'} 
B= |Bt|0);Bt e 6'} . 



The local projection of a pure state j^*) E A(E) B upon A 
is obtained by taking the partial trace over B: 



P 



(21) 



where p^- = Y^k<y-ikOi*f.- The locally projected state 
is in general a mixed state, unless pf^ is idempotent. 

It becomes harder to distinguish two states after a local 
projection since some information is lost in the partial 
trace over B. For example, the highly mixed state p, the 
pure and separable state p', and the pure entangled state 
p" given by 

. 1 



C2C1 



ctct|0)(0|c4C2 



cIct|0)(0|c3Ci- 



C3C4 



p = 



V2 72 



H0)(0| 



C2 



IC4C3 
£4 Ci + C3 



y/2 



p = \c 



t4 



3 AT 



^/2 

C2 + C3 , 



V2 



-Ci + C4- 



C2 - C3 



^/2 ' ^ V2 

become indistinguable after a local projection onto the 
orbitals 1 and 2 ({1,2}^). It is therefore in general not 
possible to draw any conclusions about the classical or 
quantum correlations in the original unprojected state 
from its local projection. 

Quantum correlations can stem from physical electron- 
electron interaction, e.g. Coulomb interaction, or be in- 
duced by a non-local projection. A non-local projection 
is qualitatively different from a local projection in that it 
projects onto a set of many-body states, and not a set of 
orbitals. The entanglement of the projected state may be 
larger than that of the original state. A trivial example 
is when the non-local projection is defined with respect 
to a single entangled state l^*), as any separable state 
with a non-zero overlap with \9) becomes entangled af- 
ter the projection. A less trivial example is given by the 
projection operator 



-It 24, 



C24.Ci^ 



:;t At 
-14. '-2t 



C2tCU- 



(22) 



The projection with P makes the pure separable spin- 



-2^ 



)|0) entangled. 



coherent state (l/2)(c[^ 

In contrast to the first example it should be noted that 
P projects onto two separable states. Moreover, the pro- 
jection itself corresponds to a simple Stern-Gerlach ex- 
periment of a spin-1 particle with a single slit to remove 
the spin ±1 components. 
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function F is: 

input: pure state l^*}, current maximal overlap 0'„^ax 

1. ^Measurement optimization: 
Diagonalize the one-particle density matrix 

p = VDV''', where pij = (^Icjcil*), such that the 
eigenvalues Yin are sorted from largest to smallest. 

2. #Detect the vacuum state: 
if (Dii = 0) return ("I'I'I') 

3. ^Perform the measurement and add another 
one-particle detector: 

for i = 1, 2, • ■ • , number of orbitals 
if (Dii > O'^ax) then 

Oma. ^ max {0'rr,ax,F CjW , 0'^ax\ ) 

end if 
done 

4. return O'^ax 
end F 

Figure 3. A pseudocode representation of the recursive depth 
first search function F [|^), O'max]- Note that the state I*!') is 
not renormalized after each measurement, and that the initial 
value of O'rnax can be set to zero. Comments are given in italic. 



Pure separable state search algorithm 



The problem of finding the maximal overlap Omax be- 
tween an N-electron pure state j^*) and any pure separa- 
ble state can be written as 



Or, 



l^")] = max 

V 



K 



(01 E 



■6uVhi|*) 



where K is the number of orbitals and V is a unitary 
matrix. A direct numerical optimization of the overlap, 
with respect to the parameters in V, becomes quickly 
very cumbersome as the number of orbitals and electrons 
increase. A practical alternative to a brute force numer- 
ical optimization is the search function F, schematically 
presented in pseudocode form in Fig. [3l The algorithm 
mimics an experiment where a sequence of one-electron 
removal measurements are performed. It is easy to visu- 
alize the search for the optimal V as an extremization of 
the intermediate intensities at each step in the measure- 
ment process. As the intensity monotonically decreases 
at each step, it is advantageous to cast the algorithm 
as a recursive depth first search function F [\'^),0'^^^], 
where 0^„a£c the currently known maximal overlap. 

Although F [|4'),0] may not return the maximal over- 
lap for a general state, it is guaranteed to be exact if 
I 'J) is separable or represents a 2-electron system. In the 
former case j^") can be represented as a single Slater de- 
terminant. Its natural orbitals are therefore either com- 
pletely occupied or empty, and the algorithm will return 



an overlap equal to one. The latter case can be proved 
by noting that the first electron removal leads to a state 
with one single electron, which is necessarily separable. 
Hence, no intensity will be lost in the second measure- 
ment, which implies that the maximal overlap is equal 
to the maximal probability of successfully removing the 
first electron. 

The search space of the algorithm can be extended 
further, e.g. by shifting the order of the electron removal 
measurements or converting electrons to holes, but such 
development is beyond the scope of this article. 



Entropic entanglement measure 

For a general many-body density operator p, the linear 
entropy 

5i[/5]=Tr[p(i-p)] 

measures how far such a state is from being pure. It has 
been known for a long time [3] that p can be represented 
by a single Slater determinant, i.e. a pure and separable 
state, if and only if the one-particle reduced density ma- 
trix pij = Tr(/3cjci) is idempotent. The linear entropy of 
the one-particle reduced density matrix 

SL[p]^TT[p{i~p)] 

measures the deviation from idempotency, which makes 
it sensitive to both classical correlations and entangle- 
ment in p. A lower bound on the contribution from 
the entanglement can be obtained by removing the 
maximal contribution from a mixed but separable en- 
semble p' = X]E=iPil^'i)(^'il' under the constrain that 
Sl[p'] = Sl[p] and N = Tr[/5] = Tr[p']. Let us consider 
a fictitious extended system, with N electrons and QN 
orbitals, for which we define 

Q 

The concavity of implies that 

Sl{(I\ < SM]. 

Evaluating the previous expression finally leads to 

SL[p']<Tr[p]SL[p] 

An analagous construction can be made for 
M = Tr[i] - N holes and QM orbitals, giving that 

SL[p']<Tiii-p]SL[p]. 

Hence, the state p is necessarily entangled if 

El[p] = Sl[p] - min(Tr[p], Tr[i - p])Sl[p] > 0. (23) 
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However, the reversed implication is not true in general. 
Any mixed state of the form 

Q 



where = 5ij and Q larger than the number of 

orbitals, gives El[p] < 0. 

It should be noted that El and the fidelity based en- 
tanglement measure 



Ep[p\ = 1 — maxTr 



'PP' 



where p' is separable, are not equivalent even for pure 
states, in the sense there exist states \'^a) and l^";,) such 
that 

Ep[\^a){\^a\]<EA\^,){\'i',\], 
i?L[|*a>(|*a|] >i?L[|*6>(|*b|]. 

For example let us set 



1 /, 



,tAt 



72 V 



-44 |o), 



1 -X 



—44 



2 ''54 



|0), 



then for 1/2 <x < 2/3 

EF[\^a){\^a\] = l>l-X = EF[\^l,){\^,\i 
EL[\-^a){\-^a\] = K (1 - X){1 + 3x) = EL[\-^b){\'fb\]- 

The difference comes from that Ep only focuses on the 
largest overlapping pure separable state |^'), while Ep 
also takes into account the entanglement between the 
states orthogonal to j^*'). 



Entanglement and spin-coherence 

Any maximally spin-polarized N-electron state with 
spin quantum number ~ s can be decomposed as 

2s N/2+S 

i*L>=E«^.n4., n 4..4.jo>. (24) 

where |0) is the vacuum state and a.y are the coefficients 
defined above. Applying the spin-coherent operator 



(25) 



to l^*! ,5) gives 



2« +Act 

" n 4.4jo)^ 



(26) 



where we have used that 



E 



(27) 



The unitary transformation 



4U ' i 1 Ai'\A (28) 



brings Eq. ((26)) back to exactly the same form as Eq. 
(|24l) . which proves that conserves the entanglement 
of 1*1,.). 

The 2s-|-l roots of unity Am = cxp[2i7rm/(2s + 1)] sat- 
isfy the important relation 

E E cxp('2^7^!^fc^)=(2. + l)^,,, 

m— — s m— — s ^ ■'^ 

(29) 

for ah -s <p,q< s. Applying Sa,„ to |'I'f,s)(^'f,3l yields 

SA„|*f,.)(*LlsL 



- E 



2.S + 1 



s /2s 



^ , ^ Ar„(s-)" \ mjmj A-(s+)" ^ 

^ I ^ ?i! 2" / 2.S + 1 I ^ n! 2" , 

m— — s \n— / \n— / 

* I** V^'*' ,1 



= E 



(2s + 1)2 " 



m.iris ,m\, — — s 

•■s lVr-l|iTiS \/iT(S iivr-1 



2s + l 



(30) 



where N;, = 27 J(J|J(2s 1) and 



) = — ^ ' , m,). (31) 



2s 
s — rrts 



Through Eq. ([30]) the thermal paramagnetic density op- 
erator 



E Priis _ \ ^ 
2s + 1 ~ ^ 



ms— — s 



ms — — s 



2s + 1 



can be written as 



E 



7n^ — — s 



2s + 1 



(32) 



(33) 



We can now introduce a set of maximally spin- 
polarized separable states l^'f use Eq. ([5T|) to de- 
fine 
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and 



' = 1^ - 



2s + 1 



E 



Finally, under the assumption that a separable density 
operator of the same form as p' gives the maximal fidelity 
'^z Pm.s'^z in Eq. ([55]). the fidelity based entanglement measure can 

2s + 1 be written as 



The entanglement preserving properties of entails that 
p' is separable as well. From Eq. (|3T|) we also have that 



\ 3,ni'^ I ^,ms / 



(34) 



which implies that the fidelity between p[-^ and p^ be- 
comes spin-independent 



F[pL,iPlr} = Tr 



F[p'sJi]. 



(35) 



Substituting Eq. ([34]) and ([35]) into F[p',p'^] gives 

■ T 2 

=Tr 



p^pWp^ 



Tr 



\ m— — £ 



Pm^^z Prn'-^z y Pr^ 



(2s + l) 



22''(2s-t- 1) 



22''(2s-|- 1) 



\ms — — s 



2s 

s — m. 



(36) 



EfIp'^] = I 



^-Ef[pJ] 
22«(2s-hl) 



E 



2s 



(37) 



